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^ ' Abstract 
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The Coulomb dissociation method to obtain the astrophysical S-factor, 5*17(0), 
. - for the ^Be(p, 7)^6 reaction at solar energies is investigated by analysing the recently 

\ measured data on the breakup reaction ^'^'^Pb(®B,''Be p)^'^^Pb at 46.5 MeV/A beam 

\ energy. Breakup cross sections corresponding to El, E2 and Ml transitions are cal- 

culated with a theory of Coulomb excitation that includes the effects of the Coulomb 
recoil as well as relativistic retardation. The interplay of nuclear and Coulomb con- 
tributions to the breakup process is studied by performing a full quantum mechani- 
cal calculation within the framework of the distorted-wave Born Approximation. In 
the kinematical regime of the present experiment, both nuclear as well as Coulomb- 
nuclear interference processes affect the pure Coulomb breakup cross sections very 
marginally. The E2 cross sections are strongly dependent on the model used to 
describe the structure of ®B. The value of 5'i7(0) is deduced with and without E2 
and Ml contributions added to the E\ cross sections and the results are discussed. 
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The radiative capture reaction ^Be(p, 7)^6 plays a crucial role in the so called "Solar 
neutrino puzzle" . The cross sections of this reaction at solar energies {E„n ~ 20 keV where 
Ecm is the center of mass relative energy of the p +'' Be system) arc directly related to 
the flux of the high energy neutrinos (emitted in the subsequent (3 decay of ^B) to which 
the ^^Cl and Kamiokande detectors are most sensitive [1]. The direct measurement of 
this reaction, at such low relative energies, is strongly hindered because of the Coulomb 
barrier. Therefore, the cross sections measured at relatively higher values of Ecm ^-re 
extrapolated to solar energies using a theoretically derived energy dependence of the 
low energy cross sections. However, absolute values of the S-factors obtained by various 
workers by following this method differ from one another considerably [2,3,4]. 

The method of Coulomb dissociation provides an alternate indirect way of determin- 
ing the cross sections of the radiative capture reaction at low relative energies. In this 
procedure it is assumed that the break-up reaction a + Z ^ {h + x) + Z proceeds entirely 
via the electromagnetic interaction; the two nuclei a and Z do not interact strongly. By 
further assumption that the electromagnetic excitation process is of first order, one can 
relate [5] directly the measured cross-sections of this reaction to those of the radiative 
capture reaction 6 -|- x — > a -|- 7. Thus, the astrophysical S-factors of the radiative capture 
processes can be determined from the study of break-up reactions under these conditions. 
A few reactions e.g. a + d -^^ Li + 7 , a + t -^'^ Li -|- 7, and -^^N + p — N + 7 have been 
studied both experimentally and theoretically using this method [see e.g. Ref. [6] for a 
recent review]. 

A first attempt has been made recently by Motobayshi et al. [8], to study the Coulomb 
dissociation of ^B into the ^Be— p low energy continuum in the field of ^ospj-^ ^j^j^ ^ 
radioactive ^B beam of 46.5 MeV/A energy [7]. Assuming a pure El excitation, the 
Monte Carlo simulation of their data predicts a 5*17(0) = 16.7 ± 3.2 eV barn, which is 
considerably lower than the value of 22.4 ±2.0 cV barn used by Bahcall and Pinsonneault 
[8] in their standard solar model (SSM) calculations. 

In this letter we perform a more rigorous analysis of the data of Ref. [7] by using a 
theory of Coulomb excitation which simultaneously includes the effects of Coulomb recoil 
and relativistic retardation by solving the general classical problem of the motion of two 
relativistic charged particles [9]. Under the kinematical condition of the present exper- 
iment, E2 transitions may be disproportionately enhanced in the Coulomb dissociation 
process. There has been quite some debate in the literature recently [10,11,12] about the 
extent of E2 contribution to the data of Ref. [7]. We would, therefore, like to reexamine 
this issue carefully, as this has important consequences for the S'i7(0) extracted from the 
data. We also include the contributions of the Ml transition which may be important for 
Ecm ill the vicinity of 0.633 MeV which corresponds to a 1+ continuum resonant state in 
^B. Furthermore, we investigate the role of the nuclear excitations. Although the cross 
sections of the pure nuclear breakup may be small as compared to those of the Coulomb 
break-up, their interference may still have some effect on the angular distributions. The 
usefulness of the Coulomb dissociation method in extracting the reliable astrophysical 
S factor from the breakup data depends on this term having negligible influence on the 
calculated break-up cross sections. We report here the result of the first full quantum 
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mechanical calculation of the Coulomb and nuclear excitations and their interference for 
this reaction, performed within the framework of the distorted wave Born approximation 
(DWBA). This is expected to highlight whether or not the nuclear effects alter appreciably 
the prediction of a pure Coulomb breakup process. 

The double differential cross- section for the Coulomb excitation of a projectile from 
its ground state to the continuum, with a definite multipolarity of order ttA is given by 
[5,6] 

(fa ^ 1 dn^x 



dildEj ^ E-y dQ 

where a^^{Ej) is the cross-section for the photodisintegration process 'j + a ^ b + x, with 
photon energy E^, and multipolarity vr = E (electric) or M (magnetic), and A = 1,2... 
(order), which is related to that of the radiative capture process a{b + x a, + 7) through 
the theorem of detailed balance. In terms of the astrophysical S-factor, we can write 

a(b + x^a + ^) = :^^exp(-27r77(£;_)), (2) 

where 77 = , with i>, Z^^ and being the relative center of mass velocity, and 

charges of the fragments b and x respectively. 

In most cases, only one or two multipolarities dominate the radiative capture as well 
as the Coulomb dissociation cross sections. In Eq. (1) rij^x^Ey) represents the number of 
equivalent (virtual) photons provided by the Coulomb field of the target to the projectile, 
which is calculated by the method discussed in Ref. [9]. S{Ecm) can be directly determined 
from the measured Coulomb dissociation cross-sections using Eqs. 1 and 2. 

The quantal treatment of the Coulomb and nuclear excitations is well known [13]. The 
expression for the differential cross sections corresponding to a transition of multipolarity 
A can be schematically represented as 



da 

dn 



oc \l drRLAkr){F^ + F^)RL,{kfr)\' (3) 



In Eq. (3), indices i and / refer to the incoming and outgoing channels. This equation also 
involves summation over the partial waves Li and Lf. F^ and F^ denote the form factors 
for the Coulomb and nuclear excitations, respectively. The expression for Fq which is 
determined entirely by the corresponding B(EX) value, is given in Refs. [13,14]. For F^, 
we take the usual collective model expression with the value of the "nuclear deformation 
parameter" being the same as that of the Coulomb one. Rl^ and Rl^. define the wave 
functions of the relative motion in the incoming and outgoing channels, respectively, which 
are obtained by solving the Schrodinger equation with appropriate optical potentials. In 
our calculations we have used the potentials given in table 3 of Ref. [15] (set D). The same 
set of potentials were used for the incoming and outgoing channels. For a Glauber-model 
calculation of the Coulomb and nuclear breakup we refer to Ref. [16]. 



2 




Figure 1: Angular distribution for ^B+^°'^Pb— i^*B*(''Be+p)+^°^Pb reactions for an El 
transition corresponding to an excitation energy of 1.137 MeV. The dashed (dotted) curve 
represent the result of pure nuclear (Coulomb) excitation calculation. Their coherent sum 
is depicted by solid curve. The scattering angle corresponds to the center-of-mass angle 
{9cm) of outgoing ^B*. 



In Fig. 1, we show the results of our full quantal calculations for the El breakup of 
on 208p|^ target at the beam energy of 46.5 MeV/A for Ecm of 1.0 MeV. We have chosen a 
higher value of Ecm where nuclear effects may be larger; the calculation reported in Ref. 
[16] is at a very low value of Ecm (100 keV). The sohd hue shows the result obtained with 
the coherent sum of nuclear and Coulomb excitations as well as their interference. The 
dashed (dotted) line depict the cross sections for pure nuclear (Coulomb) excitation. We 
have used 5000 partial waves (Lj and Lf) in the calculation of the Coulomb excitation 
cross section. We note that the pure Coulomb excitation angular distribution is a smooth 
curve approaching zero as ^ — > (the so called adiabaticity limit). This reflects the 
semiclassical nature of the process, thus the pure Coulomb excitation process should be 
amenable to the semiclassical methods [17]. The peak in this distribution corresponds 
to an angle 9min (— where a is half the distance of closest approach in a head on 

collision) below which the adiabaticity condition sets in. On the other hand, the pure 
nuclear cross sections show the typical diffraction pattern, and are at least two orders of 
magnitude smaller than the pure Coulomb one. The important point to note is that pure 
Coulomb cross sections are very similar to full calculations. Similar results are obtained 
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also for the values of Ecm of 0.8 MeV and 0.6 MeV. Thus, nuclear effects modify the total 
amplitudes very marginally in the entire kinematical regime of the data of Ref. [7] (except 
for the region very close to 6' = 0). Thus the Coulomb dissociation method can be used 
to extract reliable the astrophysical S factor for the '^Be{p, ■y)^B reaction from this data. 

In Fig. 2, we show the comparison of our calculated Coulomb dissociation double 
differential cross sections with the corresponding data of Ref. [7] as a function of the 
scattering angle 6cm of the excited (center of mass of the ^Be+p system) for three 
values of the Ecm- The calculated El, E2 and Ml cross sections are folded with an 
efficiency matrix provided to us by Naohito Iwasa (one of the authors of Ref. [7]). This 
matrix accounts for the efficiency and geometry of the detectors (including energy and 
angular spread)^. The solid lines in Fig. 2 show the calculated El cross sections obtained 
with S-factors (5*17) that provide best fit to the data (determined by minimization 
procedure). These are (17.58 ± 2.26) eV barn, (14.07 ± 2.67) eV barn and (15.59 ± 3.49) 
eV barn at Ecm= 0.6 MeV, 0.8 MeV and 1.0 MeV respectively. 

The contributions of the E2 and Ml excitations are calculated by using the radiative 
capture cross sections (T(p+'^Be— >*B+7), given by the models of Typel and Baur (TB) 
[18] and Kim, Park and Kim (KPK) [19]. We have used as input the corresponding S 
factors averaged over energy bins of experimental uncertainty in the relative energy of 
the fragments. In Fig. 2, the dashed (dashed dotted) line shows the El (with best fit 
5*17) + E2 + Ml cross sections, with E2 and Ml components calculated with TB (KPK) 
capture cross sections . It should be noted that the contribution of Ml component is 
substantial for Ecm = 0.6, while at 0.8 MeV and 1.0 MeV it is relatively quite small. 
This multipolarity was not included in the analysis presented in Ref. [10] and also in the 
calculations of the angular distributions shown in Ref. [20]. The E2 capture cross sections 
of KPK and TB are quite different from each other, while those of Ml multipolarity are 
(within 10-15%) model independent. If the KPK model is correct then E2 contributions 
to the data of Ref. [7] is not negligible. On the other hand, with the TB model the 
calculated [E1+E2+M1) cross sections are within the experimental uncertainties of the 
data. 

The best fit ''El only" Su factors, as described above, give a S'i7(0) = (15.5 ± 2.80) 
eV barn, if we use the direct extrapolation procedure adopted in Ref. [7] together with 
the TB capture model. It should be noted that in this model the El capture from both s 
and d waves to the p-wave ground state is included. In the model of Tombrello [21], which 
is the basis of extrapolation in Ref. [7] , the d wave contributions are ignored. This can lead 
to a 10-15% difference in the extrapolated S factor. To see the effect of the E2 and Ml 
components on the extracted astrophysical S factor we refitted the data of Ref. [7] with 
(i?l(best fit) -|- E2 + Ml) multiplied by a parameter ^{Ecm), which is determined by 
minimization procedure. The best fit values of the "correction factor" ^ are (0.70 ±0.17), 
(0.73 ± 0.10) and (0.75 ± 0.15) with the KPK model and (0.85 ± 0.12), (0.88 ± 0.16) and 
(0.91 ± 0.20) with the TB model, at the relative energies of 0.6 MeV, 0.8 MeV and 1.0 
MeV respectively. The values of the parameter ^{Ecm) translated into El S factors give a 

■^In Ref. [7], the detector response was taken into account by putting the theoretical cross sections as 
input to a Monte Carlo simulation program. 
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Figure 2: Comparison of experimental and theoretical Coulomb dissociation yields (cross 
section x detector efficiency) as a function of Qcm for the Ef-m values of 0.6 MeV, 0.8 MeV 
and 1.0 MeV. Sohd hnes show the calculated pure E\ Coulomb dissociation cross sections 
obtained with best fit values of S factors as discussed in the text. The dashed and 
dashed dotted curves represent the sum of E2 and Ml contributions with latter two 
components calculated with capture cross sections given in the models of TB [18] and 
KPK [19] respectively. The experimental data is taken from Ref. [7]. 
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corrected S^riO) of 11.20 ±2.02 (14.0 ±2.45) for KPK (TB) model. Thus, with the KPK 
model, the E2 and Ml contributions reduce the "El only S-factor" by more than 25%, 
while with TB model this reduction is limited to only about 10 — 12%. 

In fig. 3, the calculation for the angle integrated cross sections is compared with the 
data. The solid line shows the results for the El breakup calculated with a constant 
astrophysical S factor of 15.5 eV barn. The dashed (dashed dotted) fine shows the sum of 
the El, E2 and Ml contibutions with latter two components calculated with the capture 
cross sections of TB (KPK). It may be noted that unlike the case in Ref. [20], we have 
used efficiencies corresponding to E2 and Ml multipolarities to fold the calculated cross 
sections of these transitions. We observe that with the KPK model the cross sections up 
to Ecm — 1.2 MeV are modified by the E2 and Ml components, while with the TB model 
only the point at 600 keV is appreciably modified. 

Thus we see that E2 corrections to the data of Ref. 7 is strongly model dependent 
and it is difficult to draw any definite conclusion about their contributions. In the past 
the KPK model has been criticised by Barker [22] for their incorrect treatment of the 
resonant contribution due to Ml and E2 transitions. On the other hand Typel and 
Baur have neglected the capture from /-wave relative initial states, which could lead to 
a smaller E2 capture cross section. The calculations of several other authors also differ 
in their predictions of the E2 cross section [23]. It is therefore important to develop a 
reliable model to calculate the E2 capture cross sections for the ^Be(p, 7)^B reactions. 

Nevertheless, measurements performed at higher beam energies are likely to provide a 
kinematical regime in which Coulomb dissociation cross section are less dependent on the 
nuclear structure model, as E2 components in this regime are appreciably weaker than 
their El counterpart. In Fig. 4 we show the ratio da/dVt[{El + E2)/E1) at the beam 
energies of 46.5 MeV/A and 200 MeV/A as a function of Ocm-, with E2 cross sections 
calculated within the KPK (dotted lines) and TB (solid lines) models for the Ecm value 
of 0.6 MeV. The El angular distributions (da/dQ) decrease as 9cm increases beyond 
6m.in while those of E2 multipolarity remains fiat (in the angular region considerd in this 
figure). Thus this ratio increases with angle; the rate of the increase being determined by 
the magnitude of the E2 component. We notice that the value of this ratio is around 3.0 
(1.2) at the beam energy of 46.5 MeV/A (200 MeV/A) with E2 cross sections calculated 
in the KPK model. This implies that at the beam energies around 200 MeV and at very 
forward angles the E2 contributions calculated even in the KPK model will introduce 
relatively small modifications to pure El cross sections (they are small in the TB model 
anyway). Therefore, higher beam energies and very forward angles are expected to be 
better suited for extracting the astrophysical S-factors independent of the uncertainties in 
the nuclear strucure model of ^B. Furthermore, the higher order effects (eg. the so called 
post-acceleration) which can distort the the relative energy spectrum of ^Be + p system 
are expected to be insignificant at these energies [20]. 

In summary, we analysed the recently measured data [7] of the breakup of ^B on ^°^Pb 
target at the beam energy of 46.5 MeV/A by the Coulomb dissociation method in order to 
extract the astrophysical S factors for the radiative capture reaction ^i?e(p, 7)^-8- We used 
the first order perturbation theory of Coulomb excitation which included both the effects 
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Figure 3: Comparison of the experimental and calculated angle integrated Coulomb 
dissociation yields as a function of fragment relative energy. The sohd curve represents 
pure El cross sections calculated with a constant S factor of 15.5 eV barn. Dashed and 
dashed dotted lines show the sum of the El, E2 and Ml components with latter two 
being calculated with capture cross sections predicted by TB and KPK respectively. The 
data is from Ref. [7]. 
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Figure 4: The ratio of differential cross sections for E\ + El transitions to pure E\ 
transition (calculated with a S factor of 15.5 eV barn) as a function of Qcm for the reactions 
of Fig. 1 at the beam energies of 46.5 MeV/A (upper part) and 200 MeV/A (lower part). 
The solid (dashed) line corresponds to the case where E2 component has been calculated 
with the capture cross sections of TB (KPK). 
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of the relativistic retardation as well as Coulomb recoil simultaneously. We considered the 
breakup due to El, E2 and Ml transitions. Full quantum mechanical calculations within 
the framework of the DWBA, were performed for the pure Nuclear breakup, pure Coulomb 
breakup and their interference and it is shown that nuclear break-up contributions affect 
the total break-up amplitudes only very marginally in the kinematical regime of the data 
of Ref . [7] . Thus the assumptions of the Coulomb dissociation method are well fulfilled 
for this case. 

The £'17(0) value deduced in our analysis considering only the El cross sections is in 
agreement with the lower values reported from the extrapolation of the direct capture 
cross sections, and with that extracted from a recent indirect theoretical method [24]. 
This is significantly smaller than the value used by Bahcall and Pinsenault [8] and Turck- 
Chieze et al. [1] in their SSM calculations. The E2 break-up cross sections can reduce 
the "El only S factors", bringing the 5'i7(0) further down. However, the magnitude of 
the E2 contributions to the data of Ref. [7] depend significantly on the nuclear structure 
model used to calculate the corresponding capture cross sections and at this stage it 
is not possible to make a definite prediction. Experiments performed at higher beam 
energies are expected to be relatively less affected by the E2 component, and therefore 
they provide a better kinematical regime where the extracted S factors are likely to be 
less model dependent. 

One of the authors (RS) would like to thank Prof. R.C. Johnson and Dr. J. A. Tostevin 
for their kind hospitality at the University of Surrey and for several useful discussions on 
the subject matter of this paper. Thanks are also due to Drs. Naohito Iwasa and T. 
Motobayashi for providing us the detector efficiencies and the experimental cross sections 
of their experiment. 
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